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ABSTRACT  This  article  puts  forward  one  type  of  guidance  option 
related  to  the  maneuver  programing  of  reentering  flight  vehicles. 

The  option  in  question  makes  use  of  a  radiation  sensing  system, 
precisely  determining  the  location  of  the  starting  point  of  the 
reentering  flight  vehicle.  Use  is  also  made  of  strap-on  type  inertial 
guidance  systems  to  control  the  impact  points  of  flight  vehicles. 

KEY  TERMS  Reentry  Vehicle,  Maneuvering  Flight,  Hybrid  Guidance, 
Programmed  Guidance,  Strapdown  Inertial  Guidance,  Map  Matching 
Guidance,  Radio  Altimeter,  Application,  Digital  Simulation 

I.  INTRODUCTION 


Maneuvering  flight  is  an  important  direction  in  the  development 
of  the  technology  of  strategic  flight  vehicles.  It  is  possible  to  use 
it  in  sudden  defense,  adding  programing,  raising  precision,  and 
recovery.  Maneuvering  flight  tests  are  capable  of  giving  complete 
checks  in  all  such  areas  as  flight  vehicle  propulsion  systems, 
materials  structure,  heat  resistance  systems,  stability,  guidance 
circuits,  and  attitude  control  systems.  It  is  capable  of  carrying  out 
experiments  on  different  technological  levels.  Options  which  are 
relatively  simple  and  e^sy  to  carry  out  are  to  make  the  reentering 
flight  vehicle,  within  the  longitudinal  plane,  fly  maneuvers  according 
to  preset  programs.  In  general,  it  is  necessary  to  complete  the 
Several  steps  below: 

1.  Going  through  precise  inertial  guidance  or  celestial-inertial 
guidance  systems,  guidance  is  carried  out  on  the  flight  vehicle  in  the 
main  power  stage,  with  controls  and  corrections  to  the  free  section 
flight  orbital  plane,  such  that  the  trajectory  is  capable  of  entering 
into  the  predicted  reentry  corridor.  As  far  as  initial  flight  tests 
are  concerned,  it  is  possible  to  make  use  of  inertial  guidance  systems 
which  already  exist,  relaxing  the  requirements  for  launch  plane 
controls  and  lateral  impact  point  deviations. 

2.  Use  is  made  of  radar  altimeters  installed  aboard  the  missile 
or  other  radiation  systems,  autonomously  making  precise  determinations 
of  the  locations  of  control  points.  In  this  is  included  precise 
determinations  of  predicted  control  altitudes  and  descent  paths  at  the 


1 


altitudes  in  question.  This  requires  a  high  altitude  altimeter  with 
an  operating  range  of  from  45  Km  to  150  Km,  providing  height 
measurement  information  in  a  dispersed  way.  Perhaps  one  opts  for  the 
use  of  other-  r-adiakion  systems,  supplying  information  on  measured 
distances  corresponding  to  fixed  markers  on  the  ground. 

3.  From  control  start  points,  that  is,  setting  out  from  the  take 
off  points  for  flight  maneuvering,  the  strapdown  inectial  guidance 
option,  which  makes  use  of  inertial  measurements  to  compose  and  set  up 
the  flight  maneuvering  programs  for  reentering  flight  vehicles,  in  the 
initial  iteration  of  flight  tests,  is  capable  of  opting  for  the  use  of 
a  pure  strapdown  inertial  guidance  plan,  and,  after  the  conditions 
become  familiar,  add  on  an  imagery  matching  terminal  guidance  system. 

II.  DETERMINING  THE  SPACIAL  LOCATIONS  OF  CONTROL  START  POINTS 

The  determination  of  the  special  locations  of  control  start 
points  is  making  use  of  the  main  thrust  section  guidance  informal; ion, 
information  on  the  measurement  of  heights  and  distances  supplied  by 
the  radar  altimeters  or  other  radiation  systems  on  the  flight  vehicle 
in  the  free  section,  and,  during  the  process  of  the  flight  of  the  45 

flight  vehicle,  autonomously  making  a  fix  of  the  precise  location  of 
control  start  points.  Speaking  in  terms  of  this  method,  it  is 
possible  to  make  use  of  information  obtained  from  isochronous 
measur erne n ts  of  height  and  distances,  perhaps  making  use  of 
information  obtained  from  time  measurements  at  equal  altitudes  and 
equal  distances.  One  takes  these  types  of  information  and  carries  out 
comparisons  with  parameters  for  standard  orbits,  precisely  determining 
the  control  start  altitude  that  was  predicted  and  the  path  of  descent 
for  the  altitude  in  question. 

1.  ISOCHRONOUS  MEASUREMENTS  OF  ALTITUDE  AND  DISTANCE 

After  winding  up  the  main  thrust  section  of  the  flight  of  the 
vehicle,  one  then  enters  into  the  passive  section  of  the  flight.  At 
this  time,  the  radar  altimeters  installed  on  flight  vehicles  begin  to 
work,  obtaining  information  on  measurements  of  height  and  distance 
corresponding  with  the  surface  of  the  ground  or  measured  to  markers  on 
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the  ground.  If  the  flight  vehicle  flies  strictly  on  a  standard  orbit, 
and,  in  conjunction  with  that,  there  are  no  errors  in  the  measurements 
of  height  and  distance,  in  that  case,  the  data  for  measurements  of 
height  and  distance  will  also  be  a  series  of  ideal  values.  Actually, 
flight  vehicles  will,  more  or  less,  deviate  from  standard  orbits. 
Radiation  systems  will  also  have  measurement  noise.  At  the  current 
time,  going  through  measurement  data  in  order  co  estimate  the 
deviation  of  burnout  point  parameters  for  the  main  propulsion  stage, 
one  can  carry  on  a  step  further  and  calculate  the  index  deviation 
associated  with  flight  vehicle  reentry  maneuver  control  start  points, 
such  as  deviations  in  path  of  descent. 

In  the  case  of  flight  vehicle  motions,  one  opts  for  the  use  of 
the  model  below 

a,  t)  a=g(x,  a) 

Here,  x  is  the  vector  of  state  composed  of  the  position  vector  F  and 
the  velocity  vector  v.  w  is  the  inertial  sensing  output  information, 
u  is  the  control  force  or  interference  force  vector. 

On  the  other  hand,  the  radiation  sensing  systems  on  the  flight 
vehicle  obtain  measurement  information 

*“*<*.  t)+  e 


In  this,  z  is  the  measurment  output.  u  is  the  measurement  noise. 

Corresponding  to  standard  orbits,  equations  are  set  up  for 
state  interference  and  measurement  interference 


Here , 


A  Idx. 


b^JL\  c=irii 


l  dm  s 


dx 


(3) 


which  are  the  corresponding  partial  derivative  matrices.  v  is  the 
noise . 
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Assume  that  the  instant  the  flight  vehicle  takes  off  is  t q,  the 

instant  of  standard  burn  out  is  ^ ^ ^legible )  '  the  *nstant  °f  actual 

burn  out  is  t,  ,  •  , ,  ^  and  the  instants  of  separated  or 

k ( illegible ) 

discrete  measurements  are  t-  (t.  >,  ?  i  =  1,2,...).  When  flight 

1  IK. 

vehicles  reach  control  start  points  for  reentry  maneuvering,  the 
terminal  index  function  is  L  =  L(x,t).  After  flight  vehicles  deviate 
from  standard  states,  corresponding  to  the  instant  when  standard 
control  begins-  t  ,  there  are  deviations  of  state  <f  x"  =  cT 

3  3 

X(t  ) .  This  being  the  case,  it  gives  rise  to  an  index  of  deviation 


<5  L  =  L,6s. 


(4) 


Here,  L  is  the  determinant  vector  which  is  formed  by  the  partial 
derivatives  of  the  components  for  the  state  'x*  versus  the  function  L. 

What  must  be  paid  attention  to  in  S  xo ( illegible )  *s  tlie  ^act  ttlat 

it  is  given  rise  to  by  the  main  thrust  phase  interference  <5~w,  and  it 

r-_ 

will  be  reflected  through  0  z  in  the  free  phase.  In  conjunction  with 
this,  it  gives  rise  to  indices  of  deviation  associated  with  reentry 
control  start  points.  Going  through  analyses  of  the  considerable 
nature  of  the  systems  and  their  degrees  of  observability,  it  is 
possible  to  know  that  cT r^  is  relatively  sensitive  to  J~w,  and  that  O 
z  is  relatively  sensitive  to  o" va .  Because  of  this,  it  is  possible 
to  make  use  of  the  inertial  measurement  S  w  from  the  main  power  phase 
in  order  to  calculate  dV  ^  and  to  make  use  of  the  measure  merit  J"x 
from  the  free  phase  to  estimate  S' vo.  As  a  result  of  this,  one 
obtains  a  general  form  of  &  xq .  Going  a  step  further,  one 
calculates  <5~L,  reconsiders  time  difference  corrections,  obtaining  the 
full  deviation  associated  with  terminal  index  functions. 

From  the  theory  of  linear  systems,  on  the  basis  of  equation  (3), 
it  is  not  difficult  to  make  use  of  r  w(t)  to  solve  for  (T'x(t^)  for 
the  instant  t  .  One  separately  obtains  Jr  ( t  )  -  The  reason  for 

3  3 

this  is  that 


ta)6i(tc)  +  r)B(T)dw(T)dr 

J  to 


(5) 


When  one  pays  attention  to  <5x(tQ)  =  0  and  when  t  >  illegible) 
w(t)  -  0. 


<T 
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Because  of  this,  one  has 


6l.=di (O  =  ['* <*><*.»  r)B(r)dm(T)dr 
J  >c 


Here,  cB  ( t  t)  is  the  system  transfer  matrix.  It  has 

e(illegible) 

the  corresponding  or  analogous  components 


af  (5f(,#)  =  r)B(r)da(r)dT 

*  J  f© 


(6) 


Her  e , 


fvv 


t)  is  the  partitioned  submatrix  associated  with 


( tg  ,  t ) . 

In  another  regard,  making  use  of  xQ,  it  is  possible  to 
calculate  motion  parameters  for  any  arbitrary  instant  in  the  free 
phase 


a£(,)  =  <*>(*,  t)B(r)6a(r)dr 
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Attention  should  be  paid  to  the  fact  that,  at  this  time,  t  > 

t.  /  *  i  n  -wi  Because  of  this,  ^~w(t)  =  0.  Then,  one  has 

k ( illegible) 


<5*(/)  =<*>(*.  U< 5*.  (?) 

Making  use  of  measurements  from  equations  in  (3),  corresponding  to  the 
instant  t^ ,  one  has 


+i>i 


(8) 


Assume  that  the  ith  iteration  of  covariance  measurements  is 
R(t^).  Then,  the  minimum  variance  for  a  number  of  iterations  of 
measurement  is  estimated  to  be 
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In  order  to  solve,  from  equation  (9),  for  tile  estimated  values  of 

v  involved,  it  is  possible  to  take  the  matrix  C(t.  )  , 

e  1**1 

t  )  and  divide  it  into  two  partitioned  matrices  C  (t.,.,,  h 

e  r  1 ( illegible  V 

<"'d  c»ut(tn^ibuh  so  ■**  to  ootiin 

C(t, )$(!■•  + 


In  conjunction  with  this,  one 
Cr(ti)  jJ"  r  .  Because  of  this 


makes  o z  .  ,  . , ,  . ,  .  . 

l ( illegible ) 

,  equation  (9)  becomes 


min 


Note 


(32*  = 


Mt>)  0 
0  /?«,) 


0 

0 


It  is  then  possible  to  solve  for  the  estimated  value  of 

v  ,,,,  . ,  ,  >  ,  which  is 

e ( illegible) 


F 


iv.~(c:r-'c.t‘c:r->z" 


Or,  it  can  be  expressed  as 

Making  use  of  equation  (4),  one  obtains 


ci(#,)^-(f«)c.(#,)  J‘  'Emt.'iirvw. 


(10) 


As  far  as  the  ith  order  interference  vector  <Tw.  introduced  by  o~ 
w"  is  concerned,  it  expresses  the  i  reintegration  across  a  certain  o~oZ 
area.  Making  use  of  interference  principle  [1]  or  integral  expansion 
principle  [2] ,  it  is  possible  to  carry  out  an  integration  expansion  on 
equation  (6).  After  that  is  done,  one  obtains  <J  .. 

( t^ )  o  r .  Again,  substituting  into  equation  (10),  it 
is  possible  to  obtain  an  expression  with  the  structure  below 


i  i 


(11) 


Attention  should  be  paid  to  the  relationships  that  exist  between  the 
selection  of  A.  ^  in  equation  (11)  and  the  integration  area  or 
interval  associated  with  If,  on  the  basis  of  the  standard  time 

interval  [tQ,  tK  ]  ,  in  order  to  calculate  A-',  the  actual 
interval  of  interference  effects  is  [tQ,  t^],  at  this  time, 
corrections  should  be  carried  out  on  The  simplest  method  is 

to  introduce  a  time  correction  adjustment  quantity.  At  this  time, 
equation  (11)  becomes 


.  '  (12) 


Here,  ^  1  and  A  are  installation  coefficients  precisely 
specified  on  the  basis  of  calculating  standard  orbits.  £.  (  Jt)  is  a 
linear  function  of  time  difference.  From  the  permissible  interference 
conditions,  one  makes  precise  determinations  going  through  numerical 
simulations  or  analogues. 


2.  MEASUREMENTS  OF  TIME  AT  EQUAL  ALTITUDES 

Due  to  a  relative  similarity  in  methods  and  principles,  here,  one 
only  makes  a  simple  derivation.  First  of  all,  one  takes  system  model 
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(1)  and  changes  it  into  corresponding  equations  of  high  degree, 
introducing  H  =  S(t)  and  equation  (1).  Dividing  them  by  each  other, 
one  obtains 

After  this,  it  is  not  difficult  to  take  the  equation  above  and  write 
it  into  a  corresponding  perturbation  equation  of  higher  order  H 


(13) 


Here,  ^TT,  5  r,  and  J~v  respectively  express  time  differences,  when  one 
reaches  altitudes  that  are  the  same,  position  differences,  as  well  as 
velocity  differences.  A  is  the  time  change  coefficient  matrix  which 
is  precisely  specified  by  the  model.  Corresponding  to  equation  (13), 
one  sets  up  conjugate  equations 


(14) 


In  this,  j\T,  A.  c  ••  A  v  arg  conjugate  variables  that  correspond 
to  £*T,  r ,  and  <J"  v.  Making  use  of  the  Bliss  formula,  from  equations 
(13)  and  (14),  one  obtains 

(15) 

Appropriately  selecting  the  terminal  conditions  associated  with 
equation  ( 14 ) 


Substituting 


into 


aqua 


ti.on  (15),  one  obtains 
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6T-lr  & 


(16) 


Now,  one  has  an  appropriate  number  of  k  individual  measurement  points, 

.  ■&  f?  It 

and  one  obtains  time  deviation  information  q  ,  Q  T2 ,  ...  O  \  • 
Because  of  this,  one  has  the  relationship  below 


- iT 

rir,  x:,- 

<57, 

- 

xr,  x:. 

L<57*J 

■ 

L  Xeh  Ivk 

*•  df>.  ■* 


(17) 


Note 


5  = 


r%'i  J«i“i 

xr.xr, 

VrklW 


In  conjunction  with  this,  one  assumes  that  the  measured  times  are  all 
equally  accurate.  Because  of  this,  one  has 


jr-ur  i:>[  *■'  ]=ur  x:xb'b)-b’ 


r  <57’, 

dr, 

dT'J 


Note 


C'=cc,c,-c,)-arj:xBrfl)"5r 
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Because  of  this,  one  has 


<57=2  Ct6T' 

I  ■  I 


(13) 


As  far  as  the  coefficient  is  concerned/  it  is  precisely 
determined  by  numerical  or  digital  simulations.  Making  use  of  this 
type  of  method,  it  is  easy,  on  the  basis  of  a  certain  number  of 
measured  values  associated  with  &  T^,  to  accurately  make  precise 
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determinations  of  the  time  to  r oh  the  predetermined  altitude. 
Because  of  this,  it  is  then  also  possible  to  autonomously  determine 
the  height. 


III.  STRAP DOWN  GUIDANCE  OF  MANEUVERING  FLIGHT 

As  far  as  designing  a  lift  force  flight  vehicle  which  carries 
with  it  control  wing  surfaces  is  concerned,  the  wing  surfaces  operate 
at  a  fixed  angle,  correspondingly  producing  a  fixed  trim  angle  of 
attack.  From  this,  the  control  forces  for  maneuvering  flight  are 
produced . 

1.  ORBITAL  MODELS 

With  regard  to  the  rotating  globe,  it  is  possible  to  set  up 
orbital  equations  associated  with  the  passive  phase  of  flights: 


F  =  f> 

j  _  5 + JiojriL  +  g  q  +ae 


(19) 


Here,'?  is  the  position  vector,  v  is  the  velocity  vector  relative  to 
the  earth.  g  is  the  vector  for  the  acceleration  of  gravity  associated 
with  the  earth.  L  and  -D  are,  respectively,  the  lift  and  drag 
vectors  for  the  flight  vehicle,  a-  is  the  centrifugal  force  for  a 

_  q 

unit  mass.  ac  is  the  Kashi  (phonetic,  possiblly  Christoffel)  force 
for  a  unit  mass. 

Assume  that  u  displays  the  directional  unit  vector  for  the  axis 
or  symmetry  or  zero  lift  axis  of  the  flight  vehicle,  v/v  is  the  unit 
directional  vector  for  the  velocity  of  the  flight  vehicle.  C ®  is 
the  derivative  of  the  lift  coefficient  versus  the  angle  of  reception 
a.  p  is  the  air  density.  S  is  the  characteristic  surface  area  of  the 
flight  vehicle.  At  this  time,  the  direction  of  the  lift  force  Lq  is 
caoable  of  being  expressed  by  the  use  of  _ 
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Because  of  this  [3] 


•  •  f.  *  ’  * 


~lspc;(£>x8)xc 

it 


(20) 


In  the  derivation,  use  is  made  of  the  assumption  that  a  is  a  small 
angle,  that  is,  sin  a  a.  flaking  use  of  the  vector  product  formula, 
with  slight  additions  to  equation  (20),  changes  yield 


r0=lspc; [(!>  •  0)s-(e-o)D] 

it 


(21) 


The  magnitude  of  drag  force  D,  besides  having  the  quantity  1/2 
2 

SpC  v  ,  also  has  the  induced  drag  force  quantity  >}  1/2 
SpC^v'a  which  is  oroduced  at  the  angle  of  attack  which  is  taken 

y 

when  in  flight.  In  the  same  way,  making  use  of  a  is  done  with  the 
assumption  of  a  small  angle.  With  the  approximation  formula  1-cos  a 
ay 2,  the  induced  drag  quantity  changes  into  nq  SpC^v  (1-cos  a). 
Attention  is  paid  to  cos  a  =  (v/v)*u.  Going  through  the  necessary 
simplifications,  one  obtains 


o  ■  eAc„+2^cl)-2i7c;(*>  •  8)]» 

2 


(22) 


Assume  that  the  vector  radius  of  flight  vehicles  pointed  toward 
the  center  of  the  earth  is  R.  Then, 
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a,  =  -5  x  (5  x  /?)  =o‘R-(*  ■  R)s 
8,=  -2  5X5 


In  this  way,  one  obtains  equation  (19),  which  is  suitable  for  free 

flight  phases  and  reentry  maneuvering  flight  phases. 

If  the  range  of  flight  maneuvering  is  relatively  small,  it  is 

possible  to  ignore  the  effects  of  the  rotation  of  the  earth,  a  and 

a  need  not  be  considered.  At  this  time,  equation  (19)  changes  to 
c 

be 
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r  =  0 


(23) 


b=  g 


+ 


I  ,-P 

m 


With  regard  to  segments  of  reentry  orbits,  they  are  as  shown 
below ; 

(1)  FREE  ENTRY  PHASE.  The  altitude  is  approximately  70-90  km. 
Flight  attitude  control  and  attitude  stability  systems  guarantee  that 
flight  vehicles  take  a  zero  angle  of  attack  on  entry.  Control  wings 
for  the  next  phase  operate  to  make  good  preparations. 

(2)  CONTROL  START  PHASE 

The  altitude  is  Approximately  35-40  km.  Control  wings  maintain 
maximum  control  forces  until  lateral  overloads  reach  maximum 
permissible  values  W  =  y  max  =  30-40y. 

(3)  PULLING  EVEN  PHASE 

The  flight  vehicle  takes  the  maximum  permissible  lateral  overload 
W  y  max  to  operate  until  orbital  pulling  level  or  evening  up  and  then 
stops.  When  it  is  finished,  it  makes  the  acceleration  instruments 
aboard  the  flight  vehicle  for  the  axial  and  normal  directions  show 

detection  values — the  relative  values  for  the  normal  and  axial 
forces — as  predetermined  values,  in  order  to  guarantee  that  the 
orbital  angle  of  inclination  is  zero  or  a  certain  predicted  angle. 

(4)  LEVEL  FLIGHT  PHASE 

Making  use  of  detected  information  from  acceleration  instruments 
as  well  as  time  signal  structured  guidance  and  control  equations,  one 
carries  out  flight  maneuvering  range  calculations  for  the  path  of 
ascent,  sending  out  dive  commands  or  searching  for  character istic  land 
forms  to  carry  out  map  matching. 

(5)  DIVE  PHASE. 

After  obtaining  dive  commands,  operating  the  control  wings, 
appropriate  negative  angles  of  attack  are  realized  in  order  to  div-i. 
This  is  done  to  guarantee  the  necessary  speed  of  fall. 
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2.  GUIDANCE  EQUATIONS 

Assuming  that  one  has  installed  three  snapdown  acceleration 

instruments  on  board  the  flight  vehicle,  it  is  possible  to  detect  the 

apparent  accelerations  along  the  various  axses ,  W^,  Wyl'  ^zl‘ 

In  conjunction  with  this,  one  has  altitude  control  systems  to 

guarantee  attitude  stability.  Going  through  operations  of  the  control 

wings,  one  maintains  the  appropriate  angle  of  attack  in  order  to 

satisfy  the  required  ratio  of  forces  between  normal  and  axial 
r  4 1 

directions  .  At  this  time,  one  has 


-  —  =rtr,lcosa  —  tf^sina 

Tfl 

—  =  JJ\1sina  +  rtry1cosa  (24) 


On  the  other  hand,  we  know  that  the  ratio  between  normal  and 
axial  forces  is  related  to  angle  of  attack  a  and  Mach  number  M. 

Because  of  this,  a  is  capable  of  being  expressed  as  a  function  of  the 
K  value  for  the  ratio  of  forces  in  the  normal  and  axial  directions  and 
T  and  \T. 

c  =  a(f ,  JC)  (25) 

Using  R  to  express  the  overall  acceleration  vector  given  rise  to  by 
forces  at  work  in  the  air,  one  has,  as  a  result, 


m 


(26) 


Giving  consideration  to  the  intstallation  of  acceleration  instruments 
and  the  small  motions  of  the  missile  body  in  attitudes  approaching 
nominal  states,  on  the  acceleration  instruments  for  the  z  axis 
direction,  it  is  possible  to  have  an  output.  At  this  time,  the 
relationship  between  R  and  the  amounts  of  acceleration  detected  for 


13 


the  various  axses  is 
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Setting  up  the  conjugate  equation  corresponding  to  equation  (30),  one 
makes  use  of  the  principles  of  interference  and  obtains 


Al*+IW05O  +  i>fP. 

<  •  l 


(31) 


Here  is  the  value  selected  for  i  repeated 

integrations  or  integrals  of  overloads  on  standard  trajectories.  Z. 

* 

L  is  the  range  deviation  on  the  course  of  ascent  for  the  course 
flown  from  the  control  start  point  to  the  objective  point  relative  to 
the  course  for  standard  altitudes.  Obviously,  it  is  made  up  of  two 
sections,  that  is,  AL  =  ZL  (  +  ^  L*.  ^  L*  is  the  vertical 

range  deviation  on  the  course  of  ascent  for  the  point  in  the  main 
powered  phase  at  which  the  engines  are  turned  off.  Z*  L"^  is  the 
descending  course  deviation  at  a  predetermined  control  start  altitude. 

Making  use  of  equation  (31),  it  is  then  possible  to  construct  control 
equations.  Attention  should  be  paid,  in  equation  (25),  to  the  precise 
determination  of  K  value  parameters.  This  is  the  key  to  guaranteeing 
that  flight  vehicles  complete  their  programs  of  flight  maneuvering. 
This  will  be  realized  by  attitude  control  circuits  from  the  pitch 
channel . 


3.  DIGITAL  SIMULATION  RESULTS 

(1)  With  Regard  to  Courses  of  Descent  Associated  With  Exact 
Fixing  of  Predetermined  Altitudes. 

If  one  simply  makes  use  of  altitude  information  measured  at  a  few 
dispersed  points  in  the  orbital  descent  phase,  below  150  krn  and  above 
45  km,  it  is  possible  to  estimate  the  deviation  values  of  courses  of 
descent  as  they  correspond  to  standard  orbital  courses  of  descent.  In 
general,  they  are  90%  -  93%  of  the  deviations  of  the  actual  orbital 
courses  of  descent  relative  to  standard  orbital  courses  of  descent. 
This  is  to  say  that  the  relative  error  with  this  method  is  7%  -  10%. 

If  one  adds  measurement  data  for  an  ascent  phase  altitude  of  150  km, 
it  is  then  possible  to  raise  the  precision  of  estimated  values,  that 
is,  compensation  effects  associated  with  methods  which  are  95%  -  97%, 
or  a  relative  error  for  methods  of  less  than  5%. 
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(2)  With  Regard  to  the  Precise  Fixing  of  the  Altitudes  of 
Anticipated  Control  Start  Points. 

Calculations  clearly  demonstrate  that,  making  use  of  information 
which  is  measured  at  times  when  altitudes  a^e  equal,  in  a  phase  of  the 
descent  which  is  below  150  km  and  above  45  km,  it  is  possible  to 
precisely  determine  the  altitude  of  control  start  points.  The 
absolute  error  for  this  method  is  0.5  meter. 

(3)  With  Regard  to  Digital  Simulations  of  Guidance  Equations. 
Opting  for  the  use  of  simplified  models.  Make  V — velocity,  6 

— trajectory  angle  of  inclination,  H — height,  L — range  or  course  of 

ascent,  a — angle  of  attack,  W  ^ — the  integral  of  overload  in  the 

longitudinal  axis  direction  on  the  body  of  the  missile,  — the 

integral  of  overload  in  the  lateral  axis  direction  on  the  body  of  the 

missile,  Wx2  and  are  respectively  the  integrals  for  and 

W  gQ — the  acceleration  of  gravity  at  th?  earth's  surface, 

C  — drag  coefficient,  and  C  — lift  coefficient.  pv’S 

x  y  q=c.  x  2m 

i  o 

is  the  drag  force  for  a  unit  mass.  Y  —  C,-^^~  is  lift  f°rce 

2m 


for  a  unit  mass.  Then,  the  orbital  equations  are 


-Q-gsinO 

#  =usin0 

L-TC+-HcaxS 

#■.,  =  ( -  Qoosa  +  Y  sin  a)/g> 
W,x  *  (0  sina  +  Koo8a)/e» 


First  of  all,  as  far  as  calculating  standard  orbits  is  concerned,  it 
includes  the  first  and  second  repeated  integrals  associated  with 
overloads  in  axial  directions  and  normal  directions  after  precise 
determinations  of  control  start  times  and  the  start  of  controls. 
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After  that,  one  changes  reentry  parameters,  still  taking  the  control 
start  time  as  the  zero  point.  Going  through  adjustments  of  the 
instant  of  dive  in  order  to  guarantee  the  flight  vehicle  dropping  on 
the  objective  point,  one  simultaneously  calculates  out  the  first  and 
second  repeated  integrals  associated  with  overloads  in  axial  and 
normal  directions  from  the  control  start  point  to  the  point  of  the 
dive.  After  this  is  done,  the  data  is  formed  into  a  table.  For 
example,  the  table  of  data  for  5  orbits  from  simulation  calculations 
is : 


Off'll 

AJP  yi 

AJP  >1 

ar 

—403.2(1 

257.  US 

-5115.340 

4312.0*1 

37.90 

—  421.(52 

2(2.141 

-55(2. 946 

4791.992 

39.3* 

-429.145 

2(5.212 

-6193.303 

5013.0(9 

40.00 

129 

26(444 

|  —54(5.2(5 

521*. 174 

i 

40.53 

-451.215 

■  274.10* 

-7345.145 

57(7.614 

42.14 

One  solves  for  the  parameters 

A,t= -0.01703  A,,=  -0.001143 

1^= -0.03824  *>,  =  0. 001049 

C  =30.49427 


One  obtains  the  guidance  equation 

Mx+\,£W  yt+C=LT  ( 33 ) 

Making  use  of  (33),  one  alters  the  reentry  parameters  associated  with 
the  control  start  point.  One  simulates  a  certain  orbit.  Reentry 
speed  deviation  is  selected  as  +5  seconds.  The  deviation  of  angle  of 
reentry  is  selected  as  +0.5°,  obtaining  a  drop  point  deviation  which 
is  +500  meters.  With  precise  simulation  calculations,  solving  for 
guidance  equation  parameters,  it  is  predicted  that  it  is  possible  to 
obtain  even  better  results  in  the  control  of  drop  points. 
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